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When/(a;) is a polynomial of degree n and x^, i = 0, 1, . . ., n, are any w-f 1 points at which 
/(a;.)^0, the zeros of /(j) are known to be identical with the zeros of ^ ^ai\{x—x^, where 
ai=j{xi)l^'\xi—x^, Lucas proposed this principle for use in an electric analogue device 
for finding zeros. The present note evaluates this principle in digital computation for 
both real and complex zeros when the coefficients of /(x) are given exactly (integral or rational) 
so that the zeros of j{x) are identical with the zeros of l^Ailix — i), Ai integral. The chief 
advantages are (1) the saving of labor in tabulating i:Ai/ix-i) instead of f{x) in the neigh- 
borhood of the zero, especially for complex zeros, and (2) somewhat less work in the inverse 
interpolation for the zero. Three examples in locating a real root, and one example in 
locating a complex root were worked out in support of these findings. 



In three separate notes F. Lucas ^ describes an 
electric analogue device for calculating the roots of 
equations (also mentioned by J. S. Frame). ^ Al- 
though the principle is familiar, as far as the writer 
knows it has not been investigated from the stand- 
point of digital computation. The present note is 
intended to call attention to its advantages in finding 
the roots of polynomial equations with exact coef- 
ficients when one has a first approximation as a 
starting point. 

If f(x) denotes a polynomial of degree n, and Xt, 
i=^Oyl, . . . yn, denotes any n + l points where 
f(Xi) 9^0, from 



M = 12 ^'r'lJj\ f(^i) 



Il\Xi — Xj) 

--n{x-Xj)^ 



Jjx,) 



Ii\Xi — Xj)(x — Xt) 



there follows the well-known result (which is the 
basis of Lucas's method) that the zeros of f(x) are 
identical with the zeros of 



di 



x — Xi 



where ai= 



fix;) 



U^Xi—Xj) 



In problems where the coefficients of J(x) are 
rational, the choice of Xi=i is very convenient. By 
multiplying through, one obtains the equation in the 
form ^Ai/{x—i)=0, Ai integral, which saves a con- 
siderable number of multiplication operations, espe- 
cially when getting all or even several of the zeros of 
the same/(x). 

By choosing Xi=Xo-{-ih, for any Xq and h for which 
f(^x 1)7^0, Lucas's principle can be formulated also as 
follows: 

The zeros oi f(x), a polynomial of the nth degree, 
are identical with the zeros of A'^{J(t)/(x—t)}j m 

1 Here and elsewhere the summation is over the range to n; similarly 11 indi- 
cates a product over the same range, and n' indicates such a product with the 
vanishing factor omitted. 

2 Lucas, Compt. rend. Acad. Sci. Paris 106, 645-48 and 1072-74 (1888); lU, 
965-67(1890). 

3 J. S. Frame, MTAC 1, 337-53 (especially 347-50) (1945). 



any integer >n, where the/(^)/(x— /) is tabulated for 
any m+1 equally spaced values of t. 

Of course, in obtaining A'^{f{t)/(x—t)}, t is^the 
variable, with x as the parameter, and the resulting 
expression is then regarded as a function of x. (The 
proof is left for the reader.) 

This principle of Lucas, namely, calculating the 
zeros oi f(x) by tabulating l^ai/(;x—Xi) or XAi/{x—x)i, 
instead oi f(x) itself, in the neighborhood of a zero, 
was tried out in several difl'ereiit examples where the 
zeros had already been obtained from the tabulation 
of f{x) itself. In all examples, the use of Lucas's 
method showed a very great saving of computational 
labor. The two main advantages were (1 ) much less 
work in calculating Ai/(;x—Xi) instead of the separate 
terms of J{x), especially those of high degree,* and 
(2) somewhat less work in the inverse interpolation 
from the tabulated ^Ai/(;x—Xi) instead of the tabulated 
f{x) near the zero, which was apparent from the 
tendency of A^/A to be less in the former case.^ 

In connection with (1), in adding the separate 
terms of ^Ai/{x—Xi) consi(h^i-ably more significant 
figures were lost than in the summation of the 
separate terms oif{x); but that disadvantage is slight 
because of the ease in getting any number of places 
in Ai/(x—x^) by performing continued division on an 
ordinary 10-bank desk calculator (the x—Xi is almost 
certainly an exact number having fewer than 10 
significant figures), and the total work in tabulating 
the XAi/(x—Xi) is still much less than that in tabulating 
the/(ic). But if the x/s are not exact, or the coeffi- 
cients inj{x) are approximate, this principle of Lucas 
is severely limited in applicability, even to the point 
of not yielding a single significant figure. Thus, if in 
the example below, one were to introduce a relative 
error of 10"^^ in those exact coefficients oi f{X), it is 
apparent that the h{X) could not be obtained to even 
one significant figure. However, the choices of 
Xi=i or Xi=i—[n/2] seem most suitable for many 
problems. For example, in the case of the classical 

* There is the very well-know^n computational scheme for anX^-\-an-ix^~^-\- 
a„-22:"-2-t- . . . -f-aiar-j-ao in the form w„=a„j4-a„-i, Un~i = UnX-\-an-2, . . . , etc., 
until f(x)=Ui = U2i-\-ao, which avoids the calculation of powers of x, but which 
may be less convenient for checking and the retention of significant figures. 

5 The rate of convergence of most inverse interpolation series depends upon the 
rapidity with which the A'«/A fall off. 
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orthogonal polynomials where the coefficients are 
given exactly (for example, Hermite, Laguerre, 
Chebyshev, Legendre, etc.) the Ai/(x—Xi) can be easily 
had to any required number of places. 

As an illustration, consider the calculation of the 
zero of the polynomial 

/(X)==X6-156Z^ + 8580X4-205920Z3 

+ 2162160Z2-8648640X+8648640, 

which is near 67.28. Here/(X) was chosen to be 
His{x)l(2x), where X=4:X^ and Hi^{x) is the Hermite 
polynomial of order 13. Choosing Xi = i^ one has 



/(0)= 86 48640 
/(1)= 19 64665 
/(2) = -15 15008 
/(3)=: -27 39879 



y(4) = -24 89408 
/(5) = _13 89935 
/(6)- 69120 



The numbers Ti'{Xi—Xj) have no more than three 
digits. After obtaining the ai in their lowest terms, 
one sees that the zeros of /(X) are identical with^'the 
zeros of 



g{X)^ 



12012 392933 



94688 , 304431 



X 24(X-1) 3(Z-2) ' 4(X-3) 
155588 . 277987 . 96 



3(X-4) ' 24(X-5) ' A:-6 



which, in turn, are identical with the zeros of 
o. rv^ ^.^v^ 288288 392933 757504 

. 1826586 1244704 . 277987 2304 



X-3 



X-4. 



X-5 ' X-6 



An approximate value of the zero is X=67.2838, and 
h(X) was calculated for X=67.2838 (0.0001)67.2841. 
The separate terms oih(X) are given here to show the 
loss in significant figures upon summation: 





X=67. 2838 


X= 67. 2839 


X= 67. 2840 


X= 67. 2841 ^ 


288288/X 

-392933/ (X-l) = 

- 757504/ (X- 2) = 

1826586/ (X-3) = 

- 1244704/ (X-4) = 

277987/ (X- 5) = 

2304/(X-6) = 

h{X) = 


4284. 65693 07916 61588 7 

-5928. 03973 21819 20771 

-11603. 24613 45693 72493 6 

28414.4061178710 65494 

-19668.60397 13165 13862 9 

4463. 23120 93995 54940 5 

37. 59557 99085 56584 3 


4284. 65056 27646 43547 7 

-5928. 03078 87737 44453 8 

-11603. 22836 10507 33795 

28414. 36191 64363 08313 6 

-19668. 57289 13673 14593 4 

4463. 22404 34526 41854 5 

37. 59551 85619 71415 


4284. 64419 47565 54307 1 

-5928. 02184 53925 53255 7 

-11603. 21058 75865 44942 1 

28414.3177151390 70375 2 

-19668. 54181 15163 39043 

4463. 21687 75287 39323 1 

37. 59545 72155 86450 


4284. 63782 67673 93782 5 

-5928. 01290 20383 47054 6 

-11603. 19281 41768 05684 7 

28414. 27351 39793 51037 

-19668. 51073 17635 86746 1 

4463. 20971 16278 47235 5 

37.59539 58694 01688 2 


-0. 00000 00969 68520 


+0. 00000 00237 72288 6 


+0. 00000 01445 13214 6 


+0; 00000 02652 54257 8 



That only three values of h(X) at intervals of 0.0001 
are required in order to find the zero X to the max- 
imum attainable accuracy, follows from these 
differences: 



X 


A(X) 


A 


A2 


67. 2838 


—0.(6)09696 85200 


+0.(6)12074 08086 


+0. (12)1174 


67. 2839 


+. (6)02377 22886 


.(6)12074 09260 


.(12)1172 


67.2840 


. (6)14451 32146 


.(6)12074 10432 




67. 2841 


. (6)26525 42578 











When these differences were compared with the 
corresponding differences in j{X), it was noted 
that A7A for h{X) was only (1/25) of AVA for /(X), 
from which one can infer that the inverse interpo- 
lation is better for the function h{X). The zero 
X=67.2839 + 0.0001i? was found from the three- 
point formula^ 'p=r—r'^s, where 



Thus 



r = -2/i(67.2839)/{ A(67.2840) -A.(67.2838) } 
5 =A7{/i(67.2840)-A(67.2838) }. 



h (67.2840) -/i(67.2838) =0. (6)24148 17346, 



« H. E. Salzer, Bui. Am. Math. Soc. 50, No. 8, 513-16 (1944). 



r = -0.(7)47544 5772/0.(6)2414817346 

= -0.19688 6846, 
5=0.(12)1174/0.(6)24148=0.(6)486, 
r2=0.03876, rh=0.{^)l^M, 
-p =-0.19688 6865;X=67.28388 03113 135, 

which happens to be correct to 13 decimals. 

This method was tested upon the calculation of 
two different roots of the same tenth-degree polyno- 
mial in X (the example was to find both the smallest 
and largest zero of H2q{x)), the result being that the 
relativ^e saving of labor was ev^en greater than that 
for His(x). In fact, as the degree of the polynomial 
increases, the proportion of work saved also increases. 

Finally, this method was applied to the computa- 
tion of a complex root of a tenth-degree polynomial, 
from a rough first approximation. There the 
relative saving of labor was even greater than that 
for the examples involving real roots, due in par- 
ticular to the avoidance of the calculation of high 
powers of complex numbers. Even when the 
computation is for complex roots, the choice of 
Xi=i or Xi=i—[n/2] is still suitable, so that if the 
coefficients of f{x) are real, the calculation involves 
the sum of fractions with only real numerators. 

Washington, July 24, 1951. 
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